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Casimir energies on space-times having the fundamental domains of 
semi-regular spherical tesselations of the three-sphere as their spa- 
tial sections are computed for scalar and Maxwell fields. The spec- 
tral theory of p-forms on the fundamental domains is also developed 
and degeneracy generating functions computed. Absolute and rela- 
tive boundary conditions are encountered naturally. Some aspects of 
the heat-kernel expansion are explored. The expansion is shown to 
terminate with the constant term which is computed to be 1/2 on 
all tesselations for a coexact 1-form and shown to be so by topologi- 
cal arguments. Some practical points concerning generalised Bernoulli 
numbers are given. 
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1. Introduction. 

A recent paper [1], see also, [2], [3], on some field theory effects on half an 
Einstein Universe has prompted me to exhume and extend some previous work 
on space-times that include this case. In fact, some aspects of (free) scalar field 
theories on a hemisphere had been early discussed, [4-6], using both modes and 
images. Similar calculations were undertaken later in [7]. 

However it is not this path I wish to take here. I would rather consider the 
hemisphere as a (simple) example of a fundamental domain of a regular spherical 
tesselation and continue with the calculations begun in [8]. This reference deals 
with scalar fields, but Chang's thesis, [9], contains higher spin results and I will be 
incorporating these since the works [1,3] are concerned specifically with Maxwell 
fields. I will consider the spectral theory of p-forms in spherical tesselations and 
use it in the relevant, p = 1, case. 

Apart from its general interest, half an Einstein Universe is conformally related 
to anti de Sitter space and is of relevance to supersymmetry. 

I will concentrate on tesselations of the 3-sphere, although some results hold 
for the general sphere. The calculation is presented partly as an exercise in spectral 
geometry on bounded domains. 

I will be concerned to evaluate global, i.e. integrated, quantities such as the 
total Casimir energy, so the explicit form of the modes is not as vital as the spectrum. 

2. The geometry 

Space-time is static of the form Ix S d /T where T is a finite group of isometries 
of the sphere, in particular the complete symmetry group of a (d+ l)-dimensional 
regular polytope. The projection of this onto its circumscribing hypersphere, S d , 
yields a spherical tesselation, or honeycomb, the cells of which are the projections of 
the (i-dimensional faces of the polytope. V is generated, for example, by reflections 
in d + 1 concurrent hyperplanes. These hyperplanes (<i-flats) intersect the circum- 
scribing S d in a set of reflecting great (d— l)-spheres and divide up the honeycomb 
into a finer tesselation by forming the boundaries of |T| spherical d-simplices that 
are transitively permuted by V. One chosen such simplex can be taken as the fun- 
damental domain for this action. For d = 3 the fundamental domain is a spherical 
tetrahedron. In particular situations the tetrahedron can degenerate. An example 
is the 3-hemisphere when T has a single reflective action. 
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Although my main interest is with analysis on the fundamental domain, it is 
sometimes helpful to think of this as being the intersection of the hypersphere with 
that portion of the ambient R d+1 cut out by the reflecting hypersurfaces. It is, 
therefore, embedded in this ambient infinite alcove, hyper-kaleidoscope or Mobius 
corner, as it is variously termed, e.g. [10], [11], [12] . In the following, I sometimes 
use n = d + 1. 

The elements of V separate into two, equal sized sets distinguished by whether 
an element contains an even or an odd number of reflections. The even set is a 
(finite) subgroup of SO(d + 1) and this rotational part of V is sometimes referred 
to as the polytope group and its elements as 'direct' rotations. The whole of V 
is the extended polytope group, as it can be obtained from the rotational part by 
adjoining a reflection. This split means that the above tesselation actually consists 
of |T|/2 copies of the fundamental domain and an equal number of copies of the 
reflected fundamental domain. It is, therefore, technically not a regular tesselation. 
I term it a semi-regular tesselation. If one restricts to just the rotational part of T, 
then one does get a regular tesselation whose fundamental domain is obtained by 
sticking together a <i-simplex and its reflection. 

The 3-sphere is special because of the isomorphism SO(4) ~ SU(2) x SU(2) /Z 2 
and the rotational action of V on S 3 can be neatly effected by left and right group 
actions on SU(2) since this is isomorphic to S 3 . This structure means that the semi- 
regular tesselations of the 3-sphere are more numerous than for higher dimensions, 
e.g. [13,14]. This is, of course the same question as the number of regular polytopes, 
which goes a long way back to Schlaffi. 

3. The degeneracy problem 

In previous work, [15-17], I have considered quantum field theory on quotients 
of the 3-sphere. Since the fixed-point free action of V employed there plays no role 
in the preliminary mode analysis and in the basic construction of the (functions, I 
can, for rapidity, use the formulae in these references, with suitable modifications. 

Even with fixed points, the eigenvalues of the relevant operators are a subset 
of those on the full sphere. There are two main approaches to the calculation 
of the required degeneracies. One uses the LxR structure of V and writes the 
degeneracy as a function of the left and right 'rotation' angles which arise through 
the isomorphism, SO(3) ~ SU(2)/Z 2 , applied to each group element 7, = (7^,7^) G 
T. Since quantities, like the heat-kernel or Casimir energy, evaluated on the factored 
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sphere, i.e. on the fundamental domain, involve a group average over V (a preimage 
sum), they become just sums over the rotation angles and have to be evaluated term 
by term using explicit values for these angles. A snag is that this method only works 
nicely for the rotational part of V. Fortunately, it can be shown that the reflective 
part does not contribute to the Casimir energy for odd spheres. 

The other approach computes a generating function for the degeneracies as a 
closed form involving the degrees associated with the polytope group. This method 
is generally to be preferred but numerical agreement between the two routes provides 
a comforting check, cf [16]. 

To illustrate the techniques I look first at the scalar field and repeat some 
material from [8], which uses degrees, and then outline the angle sum approach 
which can be applied to the spin-1 case. Later, I will describe the degree approach 
for the Maxwell field after dealing with p-forms. 



4. Scalar Casimir energy by degrees 

Using the degree method allows factors of the <i-sphere to be discussed without 
many problems and I will do this as far as possible easily. 

All I say in this paper will be for the simplest case of a conformally coupled 
scalar field. Then, on the unit <i-sphere, the eigenvalues of the Helmholtz-like 
operator are well known to be perfect squares, 

f n + d -2 \ 2 

K=[ g J' n = 1 ' 3 '-'- W 

The generating function for the degeneracies, <i(/), of the eigenvalues is defined to 
be 

oo 

h(a) = J2d(l)v l , (2) 

1=0 

and of course depends only on the labelling, / = 0, 1, . . ., not on the specific form of 
the eigenvalues. 

I introduce a new parameter r by a = e~ T so that /i(a), which I also write as 
/i(t), is seen to be related to the cylinder kernel, T(r), i.e. the kernel for the square 
root of the Helmholtz-like operator, by 

T(r) = e - (d - 1)r/2 /i(r). (3) 

T(t) has a thermal significance. 
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The idea now is to construct h(r) by independent reasoning. The argument 
depends on the old fact that the eigenfunctions on the d-sphere arise from harmonic 
polynomials in an ambient Euclidean The number of these is determined by 

Molien's theorem. The harmonic generating function is, [18], the group average, 

^ (a) = ^T^det(l-aA)- (4) 

The notation is that A are the (d+1) x (d+1) matrix representatives of V considered 
as a finite subgroup of 0(d+ 1). For flexibility, a twisting, has been included, 

corresponding to the equivariant Molien theorem. x(A) = x(l) is the character in 
some representation. I will use only the simplest cases so need not enlarge on its 
general significance, see [19]. 

The usual harmonic generating function, h N (r), is when is the trivial 

character, = 1, VA Its explicit evaluation depends on Invariant Theory with the 
result that 

M-) = (i-- 2 )II 713^0' 

where the degrees d i , i=l,2,...,d+l are the degrees of the linearly independent 
generating members (basis) of the invariant polynomial ring associated with the 
action of V on the (d + l)-dimensional ambient vector space, R d+1 . Since there is 
always the basic invariant x 2 + y 2 + z 2 + . . ., of degree 2 = d d+1: say, the last factor 
on the denominator cancels the 'harmonic factor' (1 — a 2 ) leaving just the degrees 
(d 1 ,d 2 ,...,d d ) = dand 

M*)=n7i__W (5) 

i=l ^ ' 

The subscript on h N stands for 'Neumann' for the reason that all elements of 
T, including the odd ones, enter with the same sign in (4) which is the combination 
of images that gives Neumann conditions on the boundary of the physical region, 
here a fundamental domain. To obtain the Dirichlet generating function, h D (a), 
the odd elements have to come in negatively which can be achieved by setting 
equal to the sign character, det A, to get, (cf [10]), 

h ( \ 1 - a 2 ^ det (A) 

d , (6) 
11 (l_^)' 

i=l y ' 
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where d = Yli=i(^i — 1) is the degree of the Jacobian and equals the number 
of reflecting hyperplanes. I note for future reference that h N (a) is the generating 
function for a 0-form and h D (a) that for a (d+ l)-form, on IR d+1 , proportional to 
the volume form. 

Adding Neumann and Dirichlet gives the generating function when the sphere 
is factored by just the rotational subgroup of V. In this case the boundary conditions 
are periodic ones on the doubled fundamental domain. 

A very simple, but instructive, example is the relation between the full (pe- 
riodic!) sphere and the hemisphere. The generating functions for the former are 
trivially obtained by making V consist of just the identity 2 so that, in Molien's 
theorem, A is restricted to be the unit n x n matrix. Then (4) gives 

h ( \- 1 + ° 1 ° 

"sphered) (I - a )d (1 _ a )d + (1 _ a )d ' 

The two parts on the right-hand side are the Neumann and Dirichlet generating 
functions for the hemisphere for which the degrees, d, are all one, (1, 1, . . .) (see 
below) and d Q = 1. To make this plain, apply (4) to the two element group V = 
{id, R} where R is a reflection represented by the diagonal matrix 

(\ ... \ 

1 ... 
R= . . . . . 

\0 ... -lj 
A simple calculation yields the Neumann hemisphere expression if the id and R 
contributions are added and the Dirichlet one if they are subtracted. 

All the required spectral information is now in place so I just state the expres- 
sion for the total Casimir energy, E, in terms of the (^function, £(s), constructed 
from the eigenvalues (1) and degeneracies, d(l), 

E= 1 -C(-l/2). (7) 

This relation is reasonable if there are no divergences, as here. 

There is no need to extract the degeneracies from the generating function h(r) 
since the relation (3) enables the (-function to be written down immediately, 

C(s) = tF{1 2 ~ ^ [ dr (_ T) a.-i e -(d-iW2 h{r) j (8) 
J c n 



2 In this case all the degrees are one because every polynomial is invariant and the basis invariants 
can be taken simply to be the cartesian coordinates, x lt . . . ,x n , on M. n . 
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where C is the Hankel contour. 

Setting s equal to —1/2 in (8), the evaluation of E reduces to residues and 
the Taylor expansion of the integrand. Before giving the answer, I note that the 
(■-functions are examples of Barnes (-functions, [20,21], the general definition of 
which is, 

iT(l-s) f exp(-ar)(-r) s - 1 
( d (s,a\v) = — / dr 



c U i= i (! " exp(-^r)) 



2tt 

°° 

— , Re s > d . 

(a + m. u>) s 

m=0 v ' 



For simplicity, the to i are taken positive. If a is zero, the origin m = is to be 
excluded. One has the specific relations, 

C JV (a)=C d (2a,(d-l)/2|d) 
C D (a)=C d (2a,Ed i -(d-l)/2|d), 

and the calculation of the Casimir energy (7) gives generalised Bernoulli polynomi- 
als, 

g ^-w J+1 ( d+1 )! 2 n rf , <'' ((d - 1>/2|d) - 

where the upper sign is for Neumann and the lower for Dirichlet conditions. For 
odd spheres these are equal which can be seen earlier, and more generally. From 
(8), the residue gives, 

and the only contribution comes from the part of T(r) odd under r — > — r. Using 
the variable a it easy to show that the odd part is, 

a (d-i)/2(! _ a 2) *r A \ h _ (_i)ddet (A)) 

Two conclusions can be drawn from this equation. For even spheres, the sum 
reduces to one over the odd, reflective part of V. Hence the Casimir energy is zero 
on even spheres quotiented by just the rotation subgroup of V. (I only need this 
for x = 1 but it holds for all twistings.) By contrast, for odd spheres, the sum runs 
over only the rotational part of V which means that the Neumann and Dirichlet 
Casimir energies are equal (see below). This implies that on the hemisphere both 
of these are half the full sphere value since adding Neumann and Dirichlet gives the 
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full sphere quantities. The result of Kennedy and Unwin, [5], agrees with this but 
goes further because it says that the energy densities are equal. 

It is worth noting that another way of writing (12) is the inversion relation, 

T N {l/a) = {-l) d T D {a). (13) 

Such relations are used by Stanley [22]. In fact (13) is equivalent to Theorem 8.1 
in [22] and, in commutative algebra parlance, implies that the invariant polynomial 
ring is Gorenstein. 

For d = 3, (11) gives, 

E = E N = E D = B f\l | d x , d 2 , d 3 ) , (14) 

where I have used one of the basic properties of the degrees, |T| = 2 FJ d i . 

(3) 

The Bernoulli polynomial, B\ (a\d 1 , d 2 , d 3 ), and some properties, is given in 
Appendix A and leads to the following numbers for the four extended polytope 
groups, 

^3,3,3) = -J)L ^(3,3,4) 3557 



28800 ' v ' ' ' 46080 , . 

^/o „ 69391 „ /n n . 3178447 

E 3,4,3 ) = , E 3,3,5 ) = . 

v ' ' ' 414720 V ; 5184000 

It is straightforward to compute E for higher (odd) dimensions. 

The hemisphere is an example, (q = 1), of the lune, the fundamental domain 

of the group whose degrees are (q, 1, 1, . . .), [23]. (14) gives for the lune of dihedral 

angle, n/q, 

E * = 1440 9 ' (16) 
and so, for the hemisphere, E hs = 1/480, half the full sphere value, as expected. 



5. Scalar and Maxwell Casimir energies by angles 

This approach has been used by ourselves in several previous works [15-17], 
[24,25] and so I feel I can proceed rapidly and just write down the (^function derived 
in these papers from the character expression for the degeneracies after some minor 
manipulation, 

oo 

= 2jf| £ eos^eo^ E^ (^n(3cosja - cos na cos,/?) , (17) 

a, 13 n=j 
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where h(0) = 2 and h(l) = 1. The angles a and (3 are defined by, 



a = 6 R + e L , [3 = 9 R -9 L , 

where 9 L and 9 R are the rotation angles corresponding to the left and right SU(2) 
factors in the rotational part of V. As explained previously, it is sufficient in odd 
sphere dimensions, to sum over just the rotation part of Y. The sum over a and (3 
is the sum over this part. 

The identity term a = 0, (3 = formally diverges at s = —1/2 and has to be 
treated separately. For the two spins, the identity ^-functions are, 

C d (s) = ^-(( R (2s-2)-j( R (2s)), j =0,1, (18) 

which differs from the full sphere expression only by the 1/|T| volume factor. The 
non-free action of T means that other group elements have also to be treated sepa- 
rately. It turns out that it is sufficient to extract those elements with a = 0, (3 ^ 
and those with 9 L = 0, 9 R ^ 0. The first condition means that the element fixes 
a rotation plane (2-flat), the four-dimensional analogue of a rotation axis in three 
dimensions. The element belongs to the isotropy subgroup (the adjoint group) of 
r = SU(2)x SU(2), defined by the action g -> f^f -1 on g e SU(2) which group 
leaves fixed only the poles, +1 and —1, the intersections of all the fixed S^son S 3 . 

For computational completeness, the contributions of these special group ele- 
ments to the Casimir energy are written out, 



Egg = ^Qcosec 4 /3/2 - ^cosec 2 /3/2) 



E (0J) = (0£) + 1 



3=1 J=0 t 12 |p| 

E (m = ^ _ i cosec4/3/2 + 8jiC08ec 2p/<^ . 

The next step is to find, for each polytope group, the conjugacy classes with 
their orders, sizes and the values of a and (3. The group is then decomposed 
into these classes and the group sum over a and (3 performed class by class. The 
possible classes have been determined by Hurley, [26], in his computation of the 
crystal classes in four dimensions. Using his notation the class decompositions of 
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the rotation groups are, [9] , 



I ® 15E © 20K © (2 x 12) L' 
I © /' © (2 x 24)A © (2 x 6)L> © (6 + 12 + 2A)E 
®32(K®K')®12(R®R') 

I © /' © (2 x 72)A © (2 x 48)C © (2 x 6)D © (18 + 72)£ 
© (2 x 32) (K © if) © 36(R © i?') © (2 x 8) (S ® S') 
I © /' © (2 x 600)C © (2 x 30)1) © 450E © 400(K © if') 
®(2xl2)(W®W'®X® X') © (2 x 240)(Y © Y' ® Z ® Z') 
© (2 x 360)(a©6) . 

Most of this information can be found in Hurley, [26], except that for {3, 3, 5}, which 
is not crystallographic and so the new classes U, V, W, X, Y, Z and their dashed 
counterparts with opposite trace have been introduced, as well as classes a, b. 

The classes occurring in the above decompositions are given below, the notation 
being that C(x, a, det , a/iv, f3/n) refers to class 'C with invariants defined by, 

det (Al - A) = A 4 - x(A)A 3 + a(A)X 2 - X (A)X + det A , AgCcT, 

and the angles a and (3 have already been defined. Hurley does not give the values 
of a and ft, see [9]. 

7(4,6,1,0,0), J'(4, 6, -1,1,1), 

A(0, 0, 1, 1/4, 3/4), C(0, -1, 1, 1/6, 5/6), D(0, 2, 1, 1/2, 1/2), E(0, -2, 1, 0, 1), 
K(l, 0, 1, 0, 2/3) , K'(-l, 0, 1, 1, 1/3) , L(l, 1, 1, 1/5, 3/5), L'(-l, 1, 1, 4/5, 2/5), 
R(2, 2, 1, 0, 1/2), R'(-2, 2, 1, 1, 1/2), S(2, 3, 1, 1/3, 1/3), S'(-2, 3, 1, 2/3, 2/3) . 

Hurley also gives the decompositions of the extended polytope groups. For 
comparison I just exhibit that for {3 3 }, of order 120, from his Table 2.b, 

{3, 3, 3} = 1 © 15E © 24L' © 30F © 20iV © 10T' . 

The last three classes form the odd part of the group and, perforce, contribute 
nothing to the scalar Casimir energy, a fact not obvious in this approach. 

The calculation of the scalar Casimir energies via these angle forms yields the 
values obtained before, (15), as should be. 



{3,3,3} = 
{3,3,4} = 

{3,4,3} = 

{3,3,5} = 
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It turns out that the odd group part goes out for spin-one as it does for spin- 
zero. Essentially this is owing to a cancellation between E and H and will be 
detailed later. Then the Maxwell Casimir values computed again via (7), (17) and 
(18) are, 

£i(3,3,3) = ^, ^(3,3,4) 791 



14400' 1V ' ' ' 2880 

N 10453 , s 309421 ^ ' 

£,3,4,3) = , E,(3,3, 5 ) = . 

n ' ' ' 25920' n ' ' ) 324000 

Incidentally, it is possible to give a simple check of the numbers in the de- 
composition (19) which I exemplify in the {3 3 } case. Define the Shephard-Todd 
numbers by saying that b r elements of V fix a (d + 1 — r)-flat but no flat of higher 
dimensions. They are related to the degrees by Solomon's theorem, 

d+1 d+1 

n(i+(^-i)o=EM r - 

i=l r=0 

The degrees are algebraic while the Shephard-Todd numbers are geometric. 

For {3 3 }, b = 1, b 1 = 10, b 2 = 35, b 3 = 50, b 4 = 24, d x = 3, d 2 = 4, d 3 = 5 and 
these values agree with the class decomposition since the elements with a = 0, (3 ^ 
total 15 + 20 = 35 and those with a ^ 0, f3 ^ total 24 which is the number of 
those elements that fix a point and nothing else. Only even 6's are relevant for the 
rotation group. 

I just record the (6 l7 b 2 , b 3 , b 4 ) for the other cases: (16, 86, 176, 105) for {3, 3, 4}; 
(24, 190, 552, 385) for {3, 4, 3}; (60, 1138, 7140, 6061) for {3, 3, 5}. 



6. Invariant p— form theory 

As a prelude to Maxwell theory, I will look at p-forms, the eigenproblem for 
which on spheres is quite standard. One way of extending this to tessellations is 
to adapt the technique of Gallot and Meyer, [27], which is the form analogue of 
the classic M n embedding harmonic polynomial method with the coefficients of the 
ambient forms (in a local cartesian basis) being homogeneous harmonic polynomials. 
Finite 0(4) subgroup invariant theory can then be applied to these coefficients. If 
the full r is employed, the fundamental domain has a boundary and the form is 
conventionally chosen to satisfy absolute or relative boundary conditions, as will be 
clarified later. 

The extension of Molien's theorem to forms is used in the basic paper by 
Solomon, [28]. Flatto, [29], gives a summary. Although I am most concerned with 
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d = 3, in order to render the discussion a little more general, I will choose at 
the appropriate moment, for odd d, the form order p = (d — l)/2, so that the 
eigenvalues of the de Rham Laplacian for coexact p-forms on the sphere are again 
perfect squares, cf [30], [31], 

fi( P J) = (l+P + l)\ 1 = 0,1,..., (21) 

the general expression being, 

X CE (p, l)=(l + (d+ l)/2) 2 - ((d - l)/2 - p) 2 , / = 0, 1, . . . . (22) 

The eigenvalues on the factored sphere are the same, only the degeneracies change. 
A p-form, 

n 

a= ^2 ( x ) dx^ 1 A ... A dx^p , (23) 

1 

is said to be invariant under r C O(n), if for 7 e T acting on a; 1 , . . . , x n (written 

7x), 

a = 70: = a Ui ^ (7a;) dx^ 1 (7 _1 a;) A ... A dx^p (7 _1 a;) , V7 , 

or 

V... Mp (7z) = ^...^(x) ^(7) . . .^(7) , 

where ^(7) is the n x n fundamental (vector) representation of 7. 

A twisting, corresponding to (4), could be introduced by requiring a = x(l) l a - 
I will be concerned only with the trivial and the sign twistings, x{l) = 1 an d 
x{l) = det 7. In the latter case the forms are sometimes referred to as anti-invariant. 

When T is just {id} every function is an invariant function. In particular every 
homogeneous polynomial is invariant. This is the situation discussed by Gallot and 
Meyer, [27], appropriate for the full sphere embedded in M n . In the following I will 
sometimes refer to this as the 'full sphere' or 'trivial' case. 

For a general reflective group, the role of the cartesian coordinates is played 
by the polynomials, (I 1 , . . . , I n ), of the invariant basis and it is a theorem that any 
invariant p-form, (23), can be written as, 

« = Ev..,/ FlA - A ^ p ' ( 24 ) 

{M} 

where the coefficients are polynomials in I , ...,I n , written a* G k[I , . . . , I n \. 
which stands for the ring of polynomials generated by the P over the algebraic field 
k (here just the real or complex numbers) . 
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In the standard method, [27], it is shown that the form a is harmonic if the 
cartesian coefficients, a^}, are harmonic scalar functions. This is because M n is fiat 
and the differentials dx^ are parallel transported. The scaling properties of the form 
are then used to determine the radial derivatives which arise when the difference 
between the ambient and intrinsic Laplacians is computed. The label '/' in (22) is 
the degree of the polynomial coefficients in (23). The form also picks up a scaling 
factor of 'p' from the differentials in (23) (see [27] §11.9). 

The important fact is that exactly the same holds on the factored sphere, 
leading again to (22). The only difference is that (23) has to take the invariant 
structure (24). 

Because of this theorem, the differentials, 

(J 1 )** . . . (I n ) k - dTi A ... A dPv , k % e Z , (25) 

form a basis for the space of invariant homogeneous p-forms of degree /. One has 
to remember that it is the cartesian a* which are homogeneous degree /, not the 
a*. These coefficients are connected by the usual tensor relation under change of 
coordinates and so the degree, /, of the a* corresponding to the basis (25) is easily 
seen to be, 

/ = k t d t + ... + k n d n + m ix + . . . + m ip , m i = d i -\. (26) 

A basic ambient, M n , result concerns the dimension, d(p, /), of the vector space 
of (unrestricted) p-forms whose coefficients are homogeneous, polynomials of degree 
I invariant under T. The generating function is, [29] Theorem 3.16, [28], following 
from (26), 



d a (p, a) = J2 d ^ ^ = { l- a i 1 )[..(l- a < : 



e p (A 1 (A),...,A n (A)) 



(27) 



|r| ^ det(l-aA) 



where e p (x 1 , . . . , x n ) is the p-th elementary symmetric function in x 1 , . . . , x n and 
the \ are the eigenvalues of A. The m i = d i — 1 are the exponents of the group 

r. 

In view of (13), I investigate the behaviour of T a (p,a) = a n l 2 d a {p, a) under 
a — > 1/cr. A simple calculation reveals that, 

T a (p,l/a) = (-l) n T a (n-p,a), (28) 
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which is seen to be an ambient duality statement and can be re-expressed by defining 
the relative generating function of forms anti-invariant under T, 



d riP, °) = W[ E det(l-aA) = S d r(P, - > (29) 

with also 

T r (p,l/<j) = (-l) n T r (n-p,*). (30) 

Equations (28) and (30) are true generally. For reflective groups, more particularly, 
since A i = 1/Aj, the invariant and anti-invariant quantities can be related, 

T r (p, a) = (-l) n T a (p, 1/a) = T a (n - p, a) . (31) 

In the trivial case I note the trivial fact that T r = T a and the subscript can be 
omitted. 

As a special case I can check an earlier remark concerning the Dirichlet function 
by setting p equal to n in (27) to obtain, 

d a (n,a) = a m i +m *+-+ m nd a (0,a), (32) 

which is just (6) since J2i m i = ^o- 

The conclusion is that an invariant n-form corresponds to an anti-invariant 
0-form i.e. a pseudoscalar, as is well known. 

One can extend this, again to a well known statement, that, by duality an 
invariant p-form corresponds to an anti-invariant n — p-form. The boundary condi- 
tions also switch. This will be encountered later. It will turn out that invariant and 
anti-invariant forms correspond to forms obeying absolute and relative boundary 
conditions, respectively. The notation used above reflects this. 

7. p— form degeneracies 

The degeneracies on the sphere have been computed by Ikeda and Taniguchi, 
[32], and Iwasaki and Katase, [33] and I will start with this 'trivial' case. Another 
interesting discussion is provided by Week and Witsch, [34] esp. §3 and appendix, 
using ordinary tensor analysis, who derive corresponding results and compute the 
various degeneracies. The work of Paquet, [35] is a useful elaboration and part 
correction of [27]. Some other relevant discussions of p-forms on spheres and on 
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particular manifolds with boundary, such as the ball, are in Vassilevich ,[36], [37], 
Elizalde et a/, [38], [39], Copeland and Toms, [40], Camporesi and Higuchi [31], 
Cappelli and d'Appollionio, [41], from quantum field theory angles. See also Rubin 
and Ordonez, [42]. 

By restricting A to be just the unit matrix, one gets the trivial expression, 

d(M)= (T) (i-^ ' (33) 

and so (or otherwise), 

dfcO^dimP,^;')^ 1 ). (34) 

In this section I evaluate the various mode degeneracies, like (34) , on the spher- 
ical factors. In the following section, I give the equivalent generating functions, like 
(33), which are more useful to me. The reason for this doubling of effort is trans- 
parency, to make contact with other calculations and to act as a check. 

I define , the set of harmonic, homogeneous p-forms of degree /, and also 
Hf, the similar module of coclosed, harmonic forms, 

Hf = H p l nker5, 

in terms of which the required closed degeneracy on the sphere is, [33], [32], 

dim (Hf nkerd) , 

as brought up later. Here, d and 5 are the ambient derivatives. 

It can be shown, in the present polynomial case on M n , that closed implies 
exact, unless I + p = 0, and coclosed coexact, unless n — p + I = , and so there 
exists the exact sequence generated by d, 

n , p0 d , j , pp-l d pP d pp+ l d d j-, p +l „ 

(35) 

Restricting to the harmonic sub-modules in (35) gives the exact sequences, 

u — > n l+p > . . . ► n l+1 > n l > n l _ 1 > . . . ► n — > u 

(36) 

and 

Q^ H ° +p ^ ... H p-i H f 11?+? ... -1-,^p+ i _0. 

(37) 
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I treat these sequences in turn. All give recursions of the same form. Sequence 
(35) yields the relations for the unrestricted degeneracies, 

dim if = dim (Pf n ker d) + dim n ker d) , 

written, 

d b (p, I) = dg(p, l) + dg(p+l t l-l). (38) 

The recursion can be solved by either increasing or decreasing p. The first route 
gives, 

v 

d?(P, = E(- 1 )'" 1 d b(P -M + 3) + S b J P o S l0 , (39) 
i=i 

where the final term arises from the p = zero mode (a 0-form is automatically 
coclosed and, if closed, it is harmonic and constant, / = 0, grad0 = 0). It has to 
be added in by hand since the exact sequence doesn't cover this case. It exists only 
for invariant functions. In boundary condition language, invariant functions give 
Neumann conditions on the fundamental domain and anti-invariant ones, Dirichlet 
and there is no zero mode in the latter case. 

I have included a subscript, 6, to indicate the invariance, b = a, or anti- 
invariance, b = r, of the forms. Being duals, I will set *a = r and *r = a. This 
extension of the trivial situation is allowed because the exact sequences are graded 
by the periodicity type. This follows from the invariant form, (25), and its anti- 
invariant partner. A simple direct demonstration is given in section 8. 

The 'higher' part of the exact sequence leads to the equivalent formula, 

n—p 

dC(p,l) = J2(- 1 ) jd b(P + 3,l-3), (40) 

3=0 

and thence the identity, (either equate (39) and (40) or set p = in these equations), 

n 

j2(-iyd b (j,k-j) = s ba 5 k0 . (4i) 

3=0 

This can be checked for the trivial case (34), 

§<M" + ;:r)Q-(v> 
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which is zero if k > 0. However, when k = the summand on the left-hand side is 
zero unless j = when it is, unambiguously, 

o )(o) = 1 - 

This defines the right-hand side. I have used the identity, 



D D=g^G)(t-;>g^(;)(. fl :i 

with D = n + k — 1 and which follows either by recursion, or from the equation, 

(l + x) D - n =(l + x) D (l 



1 + x 

We have the end values from (39) and (40), 

dg{0,l) = 6 ba 8 lo , dg{n,l) = d b {n,l). (43) 

A general check of the identity, (41), not just for the trivial case is given later 
using the generating function version of the above which is somewhat neater. 
For the next sequence, (36), the dimensions are related by, 

dimftf = dim (Wf n ker d) + dim (H^ n ker d) , (44) 

which I write, 

h b (p,l) = hC(p,l) + hC(p+ 1,1-1). (45) 

As before, this can be solved by either increasing or decreasing p. The first 
way gives, 

hf(p, I) = Et-l)'- 1 h <>(P + 3) + Ka^io , (46) 

3 = 1 

where the last term again arises from the effect of the p = mode. 
The other route gives, 

n—p 

h c h {p, 1) = £(-1)' h b (p + 3 J- j) + 5 ba 5 p0 5 l2 . (47) 
The last term is added to give the required value 3 , 



^ If this term is not added, one obtains the spurious value h^(0, 2) = —1. Indeed, calculation 
of e.g. dimi/| on M 2 according to the formula at the bottom of p. 143 of [33] gives —1. In 
my notation this corresponds to h^ c (2,2) (see later) without the addition, which equals the 
uncorrected h^(0, 2) by duality. 
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Equating these two expressions, or setting p = 0, gives the identity, 

n 

j2(- i y h h (jj-j) = \ a (s l0 -s l2 ). (48) 

j=0 

hp is not what I want however. In the trivial case, the degeneracy, h c (p : l) : 
for closed harmonic p-forms on M n , is not the degeneracy on the sphere. For an 
ambient M n form, u, there is the collar split around the sphere, 

to = lo 1 + dr A uj 2 , 

where uj 1 and u 2 are ambient forms of orders p and p — 1 respectively. 

If to is ambiently closed, the 'tangential' part lo 1 is closed on the sphere, as 
is well known and easily checked. Projecting lo to the sphere by the inclusion, i*, 
leaves just to 1 and so the degrees of freedom in the u 2 have to be removed from h c \ 
which counts all the closed harmonic w on R". The easiest way of doing this is, 
[32], to use the harmonic coclosed module, Hf. 

Define V£ as the subspace of the space of differential forms on the sphere S d 
consisting of eigenforms associated with each eigenvalue A of the Laplacian on the 
sphere. Then there exists the isomorphism between M d+1 and S d , 

Hf nkerd< — ► V? n ker d , (49) 

where A = (/ + p)(d — p + I + 1). The dimension of the right-hand side is the 
dimension of the left and is the required closed degeneracy on the sphere. 
Turning last to the sequence, (37), the dimensions are related by, 

dim Hf = dim (Hf n ker d) + dim (Hf+* n ker d) , (50) 

or, 

h c h c {p, X) = hZ cc (p, l) + h£ cc (p +1,1-1). (51) 



frCCC - IS foe fi na i quantity I want and the recursion can be solved as before. The 
low road gives 

&? CC (P, = E(- 1 ) J " 1 ti C (P -3,1 + J) + S b J p0 8 l0 + 8 br 5 pn 8 l0 , (52) 
j'=i 
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and the high road, 

n—p 

hZ CC ( P ,l) = 52(-lYh? c (p + j,l-j). (53) 

3=0 

Duality on R n , implies 

h^ c (p,l) = hg(n-p,l), (54) 

and so, from (46), 

n—p 

hb C (p, I) = J^i- 1 ) 3 ' 1 h *b(n ~ P ~ jJ + j) + S br S n _ pfi S l0 



n—p 



(55) 



(56) 



= ^("1) J 1 h b (P J) + Mn-p,0 <*Z0 • 

3=1 

Equivalently, again using duality on M n , from (47), 

p 

hb C (p, I) = $^(-l) J K b (n -p + j,l-j) + S br 8 n _ pfi 

3=0 
P 

= X^ -1 ^ h b(P j) + Vn-p,0 • 

3=0 

As usual, setting p = in (52) and (53) gives an identity 

n 

y £(-l) i ti c (j,l-j) = 6 l0 6 ba - (57) 

3=0 

Thus the problem of finding h ccc is reduced to finding the dimension, h{p,l), 
of the harmonic p-form space on IR n . This follows easily in the trivial case since 
each form coefficient is a scalar so that the total number of harmonic forms is just 
the dimension of the p-form multiplied by the number of harmonic homogeneous 
scalar polynomials which has already been found. Therefore, in the trivial, sphere 
case, 

Hp, i) 



d+1 

P 

= d(p,l)-d(p,l-2) 

and this extends to, 



d + l\ _ fd+l-2 
I J \ 1-2 



(58) 



h b (p,l)=d b (p,l)-d b (p,l-2) 

(59) 

= K b (n-p,l) , 
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in general. This is confirmed in Appendix B. The p = 
and Besson, [10]. 

The solutions (46) and (47) become, 



case is discussed by Berard 



hZ(p, I) = V(-l)^ 1 {d b (p + j) -d b (p- j, l+j-2))+ S ba S p0 5 l0 , (60) 



and 



n—p 



h?fa = E ( d *(P + * ' " 3) ~ d b(P + 3,1-3- 2)) + 5 ba 8 p0 5 l2 . (61) 

3=0 

Everything is now in place to compute h . I am not interested in producing 
the standard expression for the degeneracy on the sphere, [33], only in numerical 
evaluation, which shows agreement in this trivial case using (58). For practical 
purposes I prefer to use generating functions. 



7. p-form generating functions 

I now construct the generating functions and start with the unrestricted de- 
generacies, d(p, I). Firstly the recursion (38), which holds for the <f s, becomes 

d b (p, a) = dg(p, a) + ad ( b J (p + 1, a) . (62) 

We can now implement the recursion from this and compare with what the 
solutions (39) and (40) give. From these two equations I find 

p oo oo 

d?(P, <r) = E d b(P + 3> l + \ a E Sp+ifi" 1 

3=1 1=0 1=0 

p oo 

= E(- 1 ) J " V_J E d b(p - h iy + ^ P o 

3=1 1=3 

*< ra) '— --^ '+6 ba 6 p0 , 

(63) 
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5>iy- v-' d b ( P -j^)-j2 -r 4 m) (p - 3-, *) 



m 

j=l x m=0 



and 



n—p oo 



(64) 



3=0 1=0 

n—p oo 

3=0 1=0 

n—p oo 

= E(- 1 ) v 'E d *(p+^o^ 

3=0 1=0 
n—p 

i=o 

where, in the step from the second to third line, the fact that d h {p, I) = for / < 
has been used. These equations follow, of course, by direct solution of the generating 
function recursion, (62). 

The identity (41) becomes, or set p = in (64), 

n 

^(-l)V^ 6 (j» = ^, (65) 

3=0 

which again can be checked directly. The general invariant situation is, 

1 n 

by a basic identity. The general anti-invariant case is, 



1 n 
(1 — cr d i) ... (1 — cr d «) V J 

(g-^O-.-^-^n) (67) 

1 ' (l-a d i)...(l-a<) ' 
= 0, 

since, apart from the trivial case, m n = 1. These results confirm (65) in general. 

The same equations hold for the various harmonic modules. For the purely 
harmonic one, Ti,f, from (46), 

P / 3~ X m \ 

^(M) = D- 1 ) j " V "?^-^)-£^i m) (P-^) { ao, 
j=l V m=0 m - 0/ (°°> 
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(the final term reflects the zero mode) and 



n—p 



3=0 

where, as before, the fact that h(p, /) = for / < has been used. 
For the harmonic coclosed module, Hf, 

h C b CC {p, a) = $>ir V- (h? c (p h C b C < m > (p - J, a) ) 

j = l V m=0 ^ 0/ 

+ °~ba°p0 + °br°pn > 

(the final term reflects the zero modes) and 

n—p 

^? CC W) = £(-l)V/i? c (p + .;>). 

The relation (55) between /i c ' c ' and ft is, 

n-p , j-i m 

h? C (p, a) = x o~i K b (n -P-J,v)-J2^j ( n ~P- 3, *) 

j=l V m=0 m - 

+ °~6r°n-p,0 

n-p / 3-1 m \ 

= E (-ir X -- J k(p + or) - E ^ *S ra) (p + * *) 

i=l V m=0 m - 0/ 

+ °~&r°n-p,0 ■> 

or, equivalently, from (56), 

p 

hZ C (p, a) = J^(-l) V /i„ 6 (n - p + j» + a 2 5 br 5 n _ P)0 
i=o 
p 

= Et -1 )^ ^ ~ J' a ) + ^A-p.O • 
J=0 

Finally, from (58), the harmonic generating function is 

oo 

h b (P, *) = E ^(P' Z ) ^ = (1 " ^) °-) • 
21 



The identity (48) becomes 



^(-l)V^ 6 (j,a)=5 ba (l-a 2 ). 

3=0 



(75) 



which, in view of (74) is the same as (65). For (74), one has therefore, in general 

l-a 2 



and 



h r (p, a) 



(l-a d i) ...(l-a<) p 
1-a 2 



(76) 



(77) 



(l-a d i)...(l-a< v 

Except in the trivial case, the numerator will cancel against the last factor in the 
denominator. Technically, it is to use these expressions that the formalism has been 
written in terms of generating functions. 

In particular for the 'hemisphere', a useful special case, the generating functions 
are, from (27) and (29), 



h h a ernisphere {p,a) = j^- 

h hemisphere^ ^ = 



(l-a) d 
1 



(I -ay 



d 

p — 1 

d 

p — 1 



+ 



(78) 



As a check, adding these gives (33). 



8. Boundary conditions 



Before proceeding with the calculation, the question of boundary conditions 
arises more acutely. This is related to the conditions imposed on the form under 
the action of V. I have considered both invariant and anti-invariant forms i.e. those 
twisted by det7. As has been mentioned, for functions (0-forms), these two types 
correspond to Neumann and Dirichlet boundary conditions. In the p-form case 
they give absolute and relative conditions, [43], [44]. This can be seen as follows. 
At a frontier hypersurface of a fundamental domain choose a boundary adapted 
coordinate system either in R n or in S d , it's not too important which. Let's select 
IR n . It is easily seen that, under the reflection in the frontier, for an invariant 
form the tangential components don't change sign, while the normal ones do, and 
conversely for the anti-invariant case. At the boundary therefore, either the normal 
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components vanish or the tangential ones do. This is the definition of absolute and 
relative and I can say that absolute = invariant, relative = anti-invariant. 

In the specific physical case of the Maxwell field, E is absolute and H is relative. 
The two contributions to spectral quantities have to be combined. Looking back to 
the relevant generating functions, (27) and (29), this removes the indirect elements 
of T, leaving, as advertised, just the rotation part. Another way of deducing this is 
via duality on R 4 as discussed in the next section. 

Incidentally, and somewhat surprisingly, anti-invariant p-form theory is not 
so straightforward as the invariant one and seems to be of somewhat specialist 
interest. The problem is to find the basis corresponding to the invariant one, (25). 
For functions, the anti-invariant polynomial algebra is obtained from the invariant 
one simply by multiplying by the Jacobian skew polynomial, J, of degree d leading 
to (6). For p-forms there exists, [45], a construction of basis anti-invariant 1-forms, 
a summary of which, for a more general situation, can be found in Solomon and 
Terao [46] §1, 6, and extensive geometric information in Orlik and Terao, [47]. See 
also [48], §4. However, since I am interested only in the generating functions, I can 
use, and have done, the isomorphism between the space of invariant p-forms and 
the space of anti-invariant (n — p)-forms (see [46] (6.6) for p = 1). 

This seems the appropriate place to confirm that the boundary conditions filter 
through the exact sequences (35), (36) and (37). 

It is sufficient to concentrate on one reflective hyperplane which, for conve- 
nience, I choose to be x n = 0. Consider the p-form a, 



Separate the coordinates 4 according to the index split, \i = (i, n) where i runs from 
1 to n — 1. Then a is invariant (upper sign) or anti-invariant (lower sign) according 
to 



It is easy to transcribe this discussion into the more elegant boundary adapted collar split, but 
I don't need such generality. 



n 




(79) 




(80) 
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Exhibit da as, 

n 

da = ^ d„ a„ „ (x) dx^° A dx^ 1 A ... A dx^p 

M() <Ml <M2 • ■ • 

n (81) 

= V A, ,, „ (aO cte M o A dx^ A ... A daA . 

Mo <Mi <M2 • ■ • 
1 

It is easy to show that (3 satisfies the relations (80) with the same signs as a. To 
spell things out: if n occurs in the set [i^ {% = to p) the coefficients are of the 
types d n a i a and d i a n Using both parts of (80) shows that each of these types 
satisfies the bottom equation of (80) with the same signs that a obeys, while, if 
the p, i does not contain n, the coefficients are typically d t a^ k and obey the top 
equation of (80) with the same sign as a. 

It is generally true that the exterior derivative preserves absolute boundary 
conditions and its dual relative ones, e.g. [43]. The present special geometry goes 
beyond this in that d and 8 preserve both sorts. 

9. Casimir energy for Maxwell by degrees 

In this section I construct specific generating functions and derive some field 
theoretic consequences in the form of Casimir energies. In the quantization of 
massless anti-symmetric tensor fields (p-forms) the demands of gauge invariance 
can be met by the well known ghosts-for-ghosts procedure, [49], [40]. On a static 
space-time, TxM, the Obukhov alternating combination of forms reduces just to 
the coexact p-form on Ai plus a harmonic combination, essentially by telescopage, 
[30] 5 . There are some question marks concerning the role of the zero modes which 
come into play at finite temperature and in the computation of the effective action. 
The zero temperature vacuum energy is unambiguous and it is this I shall look at 
here. 

The relation of the coexact and closed quantities is, 

\ c (p, l) = \ CE {p -1,0, g c (p, l) = 9 ce (p -1,0, 

on the sphere, which hold even in the presence of a boundary. 

^ The collapse of an alternating sum frequently happens. Early instances in the p-form setting 
occur in [50] and [51] and termed 'telescopage' in [52]. See also [53]. Being topological in essence, 
it appears, for example, in computations of analytic torsion, [54]. 
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As the main example I look at coexact 1-forms, p = 1, on the 3-sphere, d = 3, 
since this is equivalent to Maxwell theory, and seek to find g b E (l, a) = h b cc (2, a). 
It is sufficient to compute for just one of the symmetry types because, by duality 



on IR 4 , 



*S*(1, °) = t& CC {2, a) = h C b CC {2, a) = gg*{l, a) . (82) 



Although the computational equations can be amalgamated, they are eas- 
ily subjected to symbolic manipulation, so I leave them as they are. From (70), 
h£ cc (2, a) is seen to depend on h£ c (l, a) and h% c (0, a) which I evaluate via (73) 
giving, 

h? c (0, a) = h b (0, a) , h? c (l, a) = h b (l, a) - ah b (0, a) , 

where the h b are given by (76) and (77). Just to check, I look at the full sphere. In 
this case, 

1 + cr , , N (l + a)(4- a) 

h « (0 ' ff)= (r^j3' h a&*) = - (1 _ a) 3 > 

and, from (70), 

g^ E (l,a) = h^ c (2,a) = ^^, (83) 

(1 - a) A 

which is equivalent to the standard degeneracies, g b E (l, I) = 2(1 + + 3). 
In the general case combining the ingredients yields, 

CE (J d l+ d 2 -)- a- d 2+ d 3 -f (J d l+ d 3 — (J d l+ d 2+ d Z 

This expression is due to Chang, [9], who derived it by averaging the spin-one 
character over the group. It is one of our basic results. 

An important special case is the lune with degrees (q, 1, 1). I find, 

MM) = 7: vTTx W\ ' M 1 * ") 



and 



h 



(l-a) 2 (l-a<?) ' aV ' ' (1 - <x) 2 (l - ai) 
M °' a)= (l-a) 2 (l-a*)' KM = (1 -*)*(!- a*) 



° ( °' a)= K (1 ' a)= (l-a)2(l-a^ 



^ ( °' a)= (l-a)2(l-a«)' ^ (1 ' a)= (l-a)2(l-^) 
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Calculation produces, 

l + a q ~ 1 (2-a) 



(1-<t) 2 (1-<t«) 



as also follows from (84), as a check. 
The hemisphere case is q = 1, 

CE / \ ^ — ® / \ 

9b,hemisphere\^-i ^ _ ^3 ' (°5) 

half the full sphere value, (83). 

The coexact eigenvalues are given by (21) as /) = (/ + 2) 2 with / = 0, 1, . . . 
and so the cylinder kernel is given by, 

To repeat, T is the kernel for the positive square root of the de Rham Laplacian, 

T(r) = e - T ^ d5+5d) . (86) 

Using (84), the (^-function can be obtained in terms of Barnes functions, (9), 
as the basic equation, [9], 

3 

C CS (1, s) = C 3 (2s, Srf - dt\d) - C 3 (2s, Ed|d) , (87) 
i=i 

which holds, as said, for both symmetry types, i.e. boundary conditions. 
The Casimir energy from (7) follows as (cf (11) for scalars), [9], 



-4^p-(E 5 i 3) ^l d )- S i 3) ( l d ))- ( 88 ) 



The numbers so computed agree with those found earlier by the more clumsy 
sum-over-angles technique, (20). The lune of dihedral angle ir/q, gives, (cf (16)), 

g 4 + 5 Q 2 + 30g-3 
El = 720^ ' (89) 

and for the hemisphere E 1 = 11/240 half the full sphere value, as anticipated. 
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10. Higher spatial dimensions 



The simplifying fact in the preceding that allows the (^-function to be readily 
found is really that the eigenvalues are perfect squares. This occurs for a coexact 
p-form in the 'middle' dimension p = (d — l)/2, (22). I refer to this as 'middle' 
because the corresponding closed form has order p + 1 = n/2, assuming n is even. 
One might conveniently term these forms 'self dual' but this has other meanings. 

As an example of higher dimensions, I briefly treat a middle form in a lune 
fundamental domain. 

The harmonic generating function (76) is easily shown to be, 

e>«-M,l,...,l,q) 

where 

s(-'.u.....i..)-(-; a ) + (. + ^)(-:j) + ^: 2 2 ). 

and putting this into the mill yields finally the middle dimension g^ E {p : a) in the 
form, 

P v {<t) + o*- x Q p {<t) 

9r (P, O") 



(1 - a) 2 P (1 - a<?) ' 
where P and Q p are polynomials in a. For example, 

P 2 (a) = a 2 - 4a + 6, Q 2 (a) = -(a - 4) 

P 3 (a) = a 3 - 6a 2 + 15a - 20, Q 3 (a) = a 2 - 6a + 15 

P 4 (a) = a 4 - 8(7 3 + 28a 2 - 56a + 70, Q 4 (a) = - (a 3 - 8a 2 + 28a - 56) . 

The relation is 

A sample result is the Casimir energy for a 2-form (of either symmetry) on a 
5-dimensional lune of dihedral angle n/q, obtained by straightforward residues as, 

„ , , 5q 6 + 48q 5 + 165g 4 + 220q 3 + 75q 2 + 138q - 136 
EM = 360^ • 

The hemisphere value is E 2 (l) = 103/72. 
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11. The heat— kernel expansion 



Having the (^function on S 3 /r, it is a standard matter to compute the short- 
time heat-kernel expansion coefficients which often have an independent signifi- 
cance. In the present context, aspects of polytope and fundamental domain ge- 
ometry have been related to the scalar (p = 0) coefficients, [55]. In this section I 
compute a few of these important quantities and explore some features. 

To begin, I state my present convention for the traced heat-kernel expansion 
which is 

*W~i£C»/2* n/2 , OO) 

n=0 

with labels for the specific quantity under investigation. 

The first non-trivial example I consider is C 1 / 2 which is an integral over the 
boundary, dM, of the fundamental domain, M. The general formula is given by 
Blazic et al, ('BBG'), [56], for a smooth manifold, but the result extends unchanged 
to a singular one by dimensions. I will obtain it for a one-form on S 3 /r from the 
C-function, (87), using the generic relation on a <i-dimensional manifold, 

C(s) 7 " 7 ?7 l/2 ^ ; , s->(d-m)/2, (91) 

SV ; r((d-m)/2) s- (d-m)/2 ' V ;/ V ; 

for m = 0, 1, . . . , d— 1, d+ 1, d + 3, — The known residues of the Barnes (-function 
give, 

C?j% (d) = A ( B [ 3 \0 I d) - J>f ^ I d)j , (92) 

which computes to zero. Some technicalities are given in Appendix A. This vanish- 
ing says that there is as much Dirichlet as Neumann content in the coexact 1-form. 
The fact that the surface term is zero for the Maxwell field is an old result, [57-60], 
and follows from a cancellation between the E and H modes, where these exist. 

My next example is the constant term, C d / 2 , in the expansion which plays an 
important role, e.g. [61], and takes contributions from singularities in the manifold. 
For example, in three dimensions, C 3 / 2 , being a pure boundary quantity, involves 
an integration over the 2-dimensional frontier, dAi, of the fundamental domain, 
plus an integration over the 1-dimensional intersections ('edges') of the boundary 
pieces plus the effect of the 0-dimensional vertices of M. (which can be pictured 
as a spherical tetrahedron). Because all parts of dAi, of the various codimensions, 
are geodesically embedded in S 3 , all extrinsic curvatures vanish and one is left 
with an integration of scalar curvatures plus the vertex parts. In [62] the latter 



28 



quantities were evaluated for spin zero. Since the general expression for C 3 / 2 has 
been evaluated for a general p-form (on a smooth manifold) in [56] , it is now possible 
to repeat this calculation in the form setting. I will again pursue the calculation 
only for 1-forms on S 3 /r to give a flavour of the method. 

The total p-form ^-function is a sum of coexact p- and (p — l)-forms, 



Q ot (p, s) = Cr (P, s) + O - M , (93) 

which holds on any manifold with or without boundary 6 . 

The BBG expressions are for the general p-form but I find that, rather than 
use the combination (93), it is better to extract the coexact part, which is done by 
inverting (93) in the familiar telescopage fashion. Decreasing p gives, 

s£ e m = E(-ir q cr («,«), (94) 

q=0 



or the dual form, 



CS B (p, s) = (-If ]T (-ir « Ct^q, s) , (95) 

q=d- V 

while increasing p produces, 

<j=p+i 

In the context of Maxwell theory in a bounded domain of flat 3-space, I note 
that this has recently been used, essentially, by Bernasconi et al, [66]. 
Equating (94) and (96) gives, 

d 

Et-^H^HO, (97) 



6 The theory of forms on bounded manifolds is discussed by Duff, [63] , and Duff and Spenser, [64] . 
Probably the most complete analysis is given by Conner, [65], who introduces spaces of relative 
and absolute forms denoted by L and A4 respectively. He gives a careful treatment of the Hodge 
decomposition in the presence of a boundary. The situation is not so straightforward. There 
are (at least) two orthogonal decompositions of an arbitrary form, one 'adapted' to absolute 
conditions and one to relative. The upshot, however, is that one can proceed as in the compact 
case, and (93) is justified in general. 
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which is the archetypal telescopage identity. The earlier discussion of degeneracies 
and generating functions exhibits similar identities. 

The next standard relation I require is, in a <i-dimensional manifold, 



C(p,0) = C d/2 (p)-n (p), (98) 

where n is the number of zero modes. Labels can be attached as required. For 
example, setting s to zero in 997) yields the classic Betti number relation, 

E(-i)^ /2 = E^) 

q=0 q=0 

= X b (M), 

where x a (A / i) is the Euler number, xi.M), of the manifold and x r (.M) is the relative 
Euler number, %(.M, dM). For the factored sphere, M. = S /T, 

X (M) = 1, X (M,dM) = -l. (99) 

Simple substitution into (87) gives, for both symmetries, 



C CE (i, o) = C 3 (o, Srf - d t \d) - C 3 (o, Ed|d) 

= 3^(E B f , «i d )- B ? , ( ( 'i d )). 



(100) 



which evaluates identically to 1/2 in all cases. Some computational details are given 
in Appendix A. However, there is no need for an explicit computation, except as a 
check, because there is a topological explanation as outlined below. 

Adding the (equal) values for the two symmetries gives 1 for the fundamental 
domain of the purely rotational part of the polytope group. In this paper I simply 
present this as a mathematical fact. In particular it holds for the full sphere which I 
have discussed, [30], in connection with finite temperature effects. In fact the value 
of 1 must hold for all the doubled fundamental domains by a topological argument 
given, for the sphere, in [38] . For such a domain, there are no boundary heat-kernel 
coefficients for a general p-form by general arguments, see [43] for the smooth case. 
Hence we have, for odd d, 

C b d/2 {q) + C* d b /2 {q) = , Vq. (101) 
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I will confirm later that fixed points do not spoil this identity. 

Adding (95) and (96), and now setting p = (d — l)/2 for the middle form, 

2p+l 

tf E (p,s) + (g E ( P ,s) = - (-i) p ~ncHM + cHM)- 

q=p+l 

The two (^-functions on the left are equal and yield, when s = 0, the result, for any 
M, 

2p+l 

2(C d C /f (p) " n C E (p)) = - J2 (-V P ~ q (C b d/2 (q) + Cf /2 (q) - n b (q) ~ n* b (q)) . 

q=p+l 

Applying this to a fundamental domain, the known cohomology means that 
only one of the Betti numbers contributes (at q = 2p + 1 = d) and then (101) 
produces the final answer, 

C^ 2 ip) = \{-l) p+1 . (102) 

Curiously, the value of 1/2 accords with the conjecture in Baltes and Hilf, [58], 
that the constant term in the heat-kernel expansion is shape independent for the 
Maxwell field, and equals 1/2 (see also Baltes, [67]). It is calculated to be 1/2 for a 
cylinder of polygonal cross section. Again, the cause is a cancellation between the 
E— and H-type modes. For comparison, in the 3-ball, the value is 5/8. 

Turning now to the alternative expression for C 3 / 2) I assume, in general, that 
it can be written as a smooth part plus the contributions of the edges and vertices 
of the manifold, [68], [55], [62]. For a fundamental domain, all extrinsic curvatures 
vanish which eliminates the edge effects leaving, 

C 3 /2(p) = C 3/2 (p) + V(p), (103) 

where the smooth part, C, has been given in [56] Theorem 1.2 for the general 
p-form. 

I now set s to zero in (94) and (96) for p = 1, d = 3 and use (98). I will also 
drop the 3/2 label for ease, then 

C cs (l)=C a (l)-C a (0) + l 

= c a (i)-c a (o) + i + v a (i)-v a (o) 
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and the equivalent, 

C CE (l) = C a (2)-C a (3) 

— — (105) 
= C a (2)-C a (3) + V a (2)-V a (3). 

I have selected invariant forms, i.e. absolute boundary conditions 7 . This means 
that there is a constant zero mode for p = which is reflected in the 1 in (104) 
coming from the n in (98). The formulae in [56] allow the smooth part to be 
evaluated giving, from (104) and (105), 

C CE (l) = -^\dM\ + l + V a (l)-V a (0) 

871 (106) 
= —\dM\ + V a (2)-V a (3). 

Taking into account the computed value, C CE (1) = 1/2, adding and subtracting 
the two equations in (106) gives, 

V a (l)-V a (0) = V a (3)-V a (2) (107) 

and 

3 

T,^- 1 ) qv M = ^\dM\ + i. 

q=0 

The 0-form vertex contribution, V a (0) = V N (0), has been evaluated in [62], [69], 
and equals — 1^(0) = —V r (0) which, by duality, equals —V a (3). Hence from (107) 
V a (l) = — V a (2) and one can deduce that the vertex numbers change sign with the 
boundary conditions, 

V a (p) + V r (p) = 0. 

Adding the invariant and anti-invariant quantities corresponds to using just the ro- 
tational part of the complete polytope group. The cancellations justify our previous 
arguments leading to (102). 

I simplify (106) by giving the size of the boundary, 

I C2| 

\dM\ = 2b 1 l —± , 



which is a geometrical statement since b 1 is the number of reflecting hyperplanes. 
It can be obtained by a consideration of the coefficient C 1 ^ 2 for spin zero, [55]. 



7 Relative conditions could have been chosen. 
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For the tesselations, combining the above results, I find, 



™-™-K4±4±4zH. (108) 

This vertex contribution vanishes for the hemisphere, d 1 = d 2 = d 3 = 1, which is 
correct. 

I make the assumption that the total vertex term is the sum of the individ- 
ual vertex contributions. As in [62], in the tesselation fundamental domain there 
are four trihedral vertices, each having three dihedral angles, (it/ol^ n/a 2 , n/a 3 ), 
where (a 1 , a 2 , a s) = (<li r ? 2), (p, q, 2), (p, 2, 2), (r, 2, 2) in terms of the four tessela- 
tion types, (p, q, r) = (3, 3, 3) , (3, 3, 4) , (3, 4, 3) , (3, 3, 5) , in standard notation. So, 

V a (l) - V a (0) = I(q, r, 2) + I(p, q, 2) + I(p, 2, 2) + 7(r, 2, 2) , (109) 

where the 7's are the individual vertex contributions. They are symmetrical func- 
tions of the three angles. 

Equation (109) yields only three independent equations for the four 7's in (109). 
However one can apply the shape independent value of 1 /2 for C^/f in the case of an 
ordinary polyhedron in flat space to extract information. For example the polygonal 
cylinder (cf [58]), yields the I(p : 2, 2) value as, 



/(p,2,2) = 




if, to repeat, it is assumed that the 1/2 is distributed amongst the vertices. This 
provides the extra information needed to solve for the other contributions and simple 
elimination gives (with a check), 

/(3,3,2) = i, /(3,4,2) = A, 7(3,5,2) = A. 

Together with the known values for the 0-form vertex contributions, these numbers 
allow one to find the corresponding values for the general, i.e. unconstrained, 1- 
form. 

The Barnes ^-f unction has only a finite number of poles which cease at m = 
d — 1. The only remaining, higher coefficients in the expansion (90) can be found 
from, 

(-i) k k\c k+d/2 = a-k), fc = o,i,..., 
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which gives, from (87), 



It is shown in Appendix A that this vanishes, as must be the case (at least 
for the hemisphere). Expression (110) holds for both symmetry types and adding 
these gives the value for the periodic situation which is zero, since it is known that 
the heat-kernel expansion for coexact middle forms on odd dimensional spheres 
terminates with the C^/f term. The termination is thus generalised here to the 
tesselation fundamental domains. This has finite temperature consequences which 
I put aside for now. 

Further information on the heat-kernel coefficients on the full sphere can be 
found in [70], [31] and [38]. 

12. Conclusion 

In this paper, values have been found for the Casimir energies of scalar and 
Maxwell fields (coexact 1-forms) in orbifold factors of spheres, in particular the 
3-sphere. The manifolds are the fundamental domains of semi-regular tesselations 
and the computations were performed both by direct summation over the group by 
summing over angles and by using a degeneracy generating function to construct 
the ^-function which is given in terms of Barnes (^-functions. The results were 
specialised to the hemisphere, which is a case often looked at on its own. 

I also discussed the spectral geometry of p-forms on the fundamental domains. 
It allows, for example, the degeneracy generating function for any p to be found on 
any factor S d /r. 

Aspects of the short-time expansion of the heat-kernel were discussed. In 
particular it was shown, as expected, that the surface term, C 1 / 2 , vanished. The 
constant term was also computed to be 1/2 for the coexact 1-form on all factors 
and shown to be so topologically. The expression was analysed in some detail, 
bringing out the contribution of the vertices. It was also shown that the expansion 
terminates with the constant term on the fundamental domains. 

Another path to vertex contributions is via the heat-kernel on the generalised 
cone, [71], which is a bounded Mobius corner, having as two-dimensional base a 
fundamental domain of, this time, S 2 /r, the familiar tiling of the two-sphere, [72]. 
I defer this to another time. 
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Having the (^-functions and generating functions allows other spectral quanti- 
ties to be evaluated, for example the functional determinants and also finite tem- 
perature quantities. 

Although I have looked mainly at 1-forms on factors of the three sphere, it is 
possible to analyse forms of the middle order in any dimension with very similar 
results. The (^-functions are again of Barnes type. Even for forms of any order, with 
non perfect square eigenvalues, the calculation can be pursued, with determination, 
e.g. [40], [73], [38], [74] . 

I have only computed integrated quantities. It is possible to extend the same 
techniques to the local objects such as the Green function. 



Appendix A 

I give some technical niceties connected with the generalised Bernoulli poly- 
nomials in terms of which the field theory quantities have been expressed. They 
satisfy several important identities, [75], often leading to simplifications, or at least 
variety, in their evaluation. One such is the recursion 

B^(d % | d) = fl(»>(0 | d) + d^B^iO | dj , (111) 

where d^ stands for the set of n-degrees, omitting d^ 

d; = (d l ,...,d i _ l ,d i+1 ,...d n ). 

The Bernoulli polynomials of argument, B^\o\d), are the generalised 
Bernoulli numbers, B^\d), and, as an example, the coexact 1-form vacuum energy 
(88) can be re-expressed purely in terms of numbers as 

El = " 48(^3 ( 2S * 3)(dl ' ^ ds) + Ad ^^ rf 3 )+4^ 2) K, d 3 ) 

+ 4d 3 Bi 2 \d 1 ,d 2 )). 



This is not strictly necessary as the general expression for B^ (x | d) could be found. 
In a similar fashion, the heat-kernel coefficient (92) becomes 

C V2( d ) = T ^-r(^\d 11 d 2l d 3 )+d 1 B^(d 2l d 3 )+d 2 Bi ) 2 \d 11 d 3 ) 

d!d 2 d 3 \ (n2) 

+ d 3 B ( i\d 1 ,d 2 ) ). 
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and for (100), 



^ 3 C /f (d) = ^^(2Bi 3 \d 1 ,d 2 ,d 3 ) + Sd 1 Bi 2 \d 2 ,d 3 ) + Sd 2 Bi 2 \d 11 d 3 ) 

+ 3d 3 B?\d 1 ,d 2 )^J (US) 

1 

~ 2 ' 

The expression, (110), for the higher coefficients can be re-expressed, 

= STSL ( <3 + 2fc) E-^^A) + 2B« 3 + » 2t (d)) . (114) 

The spin zero vacuum energy (11) can also be subjected to some (inessential) 
manipulation. On factors of the two-sphere, the form was given in terms of the two 
degrees, [8], 

which can also be written in terms of the Shephard-Todd numbers, b l and b 2 , using, 
[55], 

Bf ((d - l)/2 | d) = i 6 X - 6 2 + ^) (115) 

as 

= ±-jTf h - b i~l)> l r l = 2d ^ ■ 

The advantage of using the Shephard-Todd numbers is that, for a given is, 
Bl d \x\d) can be displayed as an explicit function of the dimension, d. One path 
to this representation is via Todd Polynomials, T k , as Hirzebruch, [76], gives the 
relation, 

where the c s are the elementary symmetric functions of the degrees d i (i = 1, . . . , d). 
The relations needed to turn the c s into elementary symmetric functions, the b r , of 
the exponents, m i (i = 1, . . . , d + 1) are detailed in [55]. For calculational rapidity 
here, I just use the c s and note the first three Todd polynomials, 

T\ = - C\ , T 2 = — {c\ + c 2 ) , T 3 = — c x c 2 , 
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leading to, 

B[ d \x\d)=x-^ l d l 
Bi 2 \d 1 ,d 2 ) = ^((d 1 + d 2 ) 2 + d 1 d 2 ) 
B 3 '(d 1 , d 2 , d 3 ) = — -((i-L + d 2 + d 3 ){d x d 2 + d 2 d 3 + d 2 d 3 ) , 

employed in computing (112) and (113). 

Alternatively, one can use the general formula in terms of ordinary Bernoulli 
numbers, 

Bl n) (d) = 7T~~T7 # • • • d " B ; ■ ■ ■ B s n , (H6) 

*1- • • • *n- 

where the sum is over all s x , . . . , s n , positive or zero, satisfying s 1 + s 2 + . . . + s n = v. 

In the expression for the higher coefficients, (114), the lower index on the S's is 
now bigger than the upper one and (116) is called for. Writing out B\, where v is 
odd and greater than 3, the condition s 1 + s 2 + s 3 = v means that one, and only one, 
of s 1 ,s 2 ,s 3 must equal 1 because only B x = —1/2 of the odd Bernoulli numbers 
is non-zero. Breaking up the summation by taking this into account effectively 
reduces n and v by one at the expense of producing three sums. The outstanding 
factors then lead to the identity, 

2S( 3 )(d) + ^rf l S^ 1 (dJ=0, z,odd>3, 

i 

which was used in the text to show the termination of the coexact heat-kernel 
expansion from (114). 

The same arguments can be applied for v = 3. Then there is the additional 
possibility that all three s 1 ,s 2 ,s 3 equal one. This gives the extra term, 3d 1 d 2 d 3 , 
and is another way of deriving (113). Higher n can be treated similarly. 

Finally I give a polynomial which I have used, 

Bf\x | d) = - ^- [df + d% + d\ - hd\d\ - hd\d\ - hd\d\ 

— Vad\d 2 d 3 — 150^2^3 ~~ ^>d x d 2 d\ v \ 

— x \d\d 2 + d x d\ + d\d 3 + d 2 d 3 + d x d\ + d 2 d 3 + Sd^^^ 
+ x 2 [d\ + M x d 2 + d 2 2 + 3d t d 3 + M 2 d 3 + d 2 3 ] 

- x 3 [d t + d 2 + d 3 ] + x 4 . 
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Appendix B 



I have found it constructive to check a specific numerical example. I choose 
the following parameters p = 1, / = 2 in M 3 , i.e. d = 2 (n = 3). The general 
homogeneous polynomial of second degree is the conic, 

ax 2 + 2hxy + 2gxz + 2fyz + by 2 + cz 2 , 

and is therefore 6 dimensional on the full sphere. The harmonic condition imposes 
one constraint, giving 5-dimensions. A basis would be, xy, xz, yz, x 2 — y 2 , y 2 — z 2 . 
One now has to impose the closed condition. In vector language for a 1-form, = A, 
dA = curl A and so we get 3 identities, each of which is linear in x,y,z, giving 9 
conditions not all of which are independent in view of divcurl = 0, which is one 
scalar identity, making actually 8. This leaves, finally, a 15 — 8 = 7 dimensional 
module. 

The standard expression for the closed degeneracy on S 2 derived in the cited 
references gives the value 7 for h (1, 2). 

Applying now, (53), (52), (56) and (59) in detail, I find, 

h ccc (l, 2) = h cc (l, 2) - h cc (2, 1) + h cc (3, 0) 

= (Ml, 2) - h(0, 1)) - (h(2, 1) - h(l, 0)) + (h(3, 0) - h(2, -1)) 
= (d(l, 2) - d(l, 0) - d(0, 1)) - (d(2, 1) - d(l, 0)) + d(3, 0) 



and 



= 3.6-3.1-3.3 + 1.1 = 7 

h ccc (l,2) = h cc (0,3) = h(0,3) 
= d(0,3) -d(0, 1) 

MlHl 

= 10-3 = 7. 

In the coclosed case one has 5 A = div A = which is just one condition 
linear in x, y, z, equivalent to 3 conditions so that there are 15 — 3 = 12 independent 
components and this is born out by the calculation of h c (2, 2) = 12 in 3 dimensions. 
This just reflects the duality on M 3 between a closed p-form and a coclosed (n—p) = 
(3 — p) one. 
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I can also check (59) for the harmonic degeneracies in the invariant /anti- 
invariant cases. For the hemisphere there is just one reflection, say z — > —z. The 
polynomial basis splits into the 4-dim invariant set, x 2 ,y 2 ,z 2 ,xy, and the 2-dim 
anti-invariant one, xz, yz. 

In cartesian coordinates the invariance/anti-invariance of the 1-form (vector) 
under the reflection implies, 

A x (x,y,-z) = ±A x (x,y,z) 
A y (x,y,-z) = ±A y (x,y,z) 
A z (x,y,-z) = TA z (x,y,z) , 

upper sign for invariance. Therefore, for the invariant case, A x and A y need the 
invariant polynomial basis and A z the anti-invariant one, and conversely for an anti- 
invariant form. This makes explicit the nature of absolute and relative conditions 
the tangential components on the hemisphere being linear combinations of A x and 
A y . In another language, invariant vectors are polar vectors while anti-invariant 
ones are termed axial. 
Therefore, 

d a (l, 2) = 4.2 + 1.2 = 10 , d r (l, 2) = 2.2 + 1.4 = 8 . 

For the / = case, the coefficients are constants and clearly 

d a (l,0) = 1.2 = 2, rf r (l,0) = 1.1 = 1, 

giving the combinations, 

d a (l, 2) - d a (l, 0) = 10 - 2 = 8 , d r (l, 2) - d r (l, 0) = 8 - 1 = 7 

(adding to give 15). 

These values agree with the hemisphere Poincare series, 

p a (a) = 2 + 5a + 10a 2 + . . . p r {a) = 1 + 4a + 8a 2 + . . . . 

The number of harmonic invariant/anti-invariant forms can also be counted 
directly from the 5.3 harmonic basis. Thus, by inspection, the invariant basis is 
xy, x 2 — y 2 , y 2 — z 2 and the anti-invariant one xz, yz, so 

h a (l, 2) = 2.3 + 1.2 = 8 , h r (l, 2) = 2.2 + 1.3 = 7 
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in agreement with the above. 

The dimension can be increased without difficulty to, say, 4, just to make sure. 
The full basis is 10-dim, 

x 2 , y 2 , z 2 , w 2 , xy, xz, xw, yz, yw, zw 

which splits into a (7-dim) and r (3-dim) parts, 

x 2 ,y 2 ,z 2 ,w 2 ,xy,xw,yw and xz,yz,zw 

hence 

d a (l,2)-d a (l,0) = (7.3 + 1.3)-1.3 = 21 
d r (l, 2) - d r (l, 0) = (3.3 + 1.7) - 1.1 = 15 
The total harmonic basis is 9-dim, 

x 2 — y 2 , y 2 — z 2 ,z 2 — w 2 ,xy, xz, xw, yz, yw, zw 

and this splits as 6+3, 

x 2 — y 2 ,y 2 — z 2 , z 2 — w 2 ,xy, xw,yw and xz,yz,zw 

giving 

h a (l, 2) = 6.3 + 1.3 = 21 , h r (l, 2) = 3.3 + 1.6 = 15 
again in agreement. 

It is also advisable to check another p, say p = 2. The form is 6 dimensional 
(in 4 dimensions) 3 components being symmetric and 3 anti-symmetric under the 
reflection. 

d a (2, 2) - d a (2, 0) = (3.7 + 3.3) - 1.3 = 27 

h a {2,2) = 6.3 + 3.3 = 27 
d r (2, 2) - d r (2, 0) = (3.3 + 3.7) - 1.3 = 27 

h r (2,2) = 6.3 + 3.3 = 27. 
This result also exhibits the duality in (59). 
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